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Abstract 

A unified description of the free-space cylindrical vector beams is presented, which is an integral 
transformation solution to the vector Helmholtz equation and the transversality condition. The 
amplitude 2-form of the angular spectrum involved in this solution can be arbitrarily chosen. When 
one of the two elements is zero, we arrive at either transverse-electric or transverse-magnetic beam 
mode. In the paraxial condition, this solution not only includes the known J\ Bessel-Gaussian 
vector beam and the axisymmetric Laguerre-Gaussian vector beam that were obtained by solving 
the paraxial wave equations, but also predicts two new kinds of vector beam, called the modified- 
Bessel-Gaussian vector beam. 
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A free-space cylindrical vector beam is of spatially inhomogeneous polarization that is 
rotationally symmetric with respect to the propagation axis. Due to this unique polarization 
characteristic, the cylindrical vector beam has attracted much attention in both optical 
physics flflfla and applied optics flflQ. The mathenrafeal expression for a paraxial 
indrical vector beam is usually obtained by solving the paraxial wave equation |8|, [9|, HO, 
nd the known kinds of cylindrical vector beam include the Bessel-Gaussian 
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lCj and Laguerre-Gaussian beams [12j. In this paper, we present a unified description of 



the cylindrical vector beam, from which we obtain for the first time two new kinds of vector 
beam, called the modified-Bessel-Gaussian vector beam. This beam description represents 
an integral transformation solution to the vector Helmholtz equation and the transversality 
condition. 

Consider a light beam that propagates in the positive x direction in source-free space. 
The electric-field vector E of the beam satisfies the vector Helmholtz equation, 



V 2 E(x) + fc 2 E(x) 



0. 



(1) 



subject to the transversality condition 



V ■ E(x) = 0. 



(2) 



For a rotationally symmetric beam with respect to its propagation axis, it is convenient 
to make use of the cylindrical coordinate system, in which x = xe x + r, where r = re r = 
e y r cos (f) + e z r sin <fi is the polar coordinate. It has been shown that the electric- field vector of 
;he beam can be represented by the following integral over the plane- wave angular spectrum 

E(x) = -W / A(k p , <p) exp(.k • x)dft, (3) 

^ J Jk p <k 



where k = k x e x + k p is the wavevector of the plane wave, k p = k p e p = e y k p cos (p + e z k p sin ip, 

k = (k 2 - kiyi\ 

A{k p ,ip) = PA{k p ,ip) (4) 



is the amplitude vector of the angular spectrum, 



A(k p , if) 



A, 
A 



(5) 
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H 



is the amplitude two-form [14J of the angular spectrum 
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(6) 



is a 3 x 2 matrix that plays the role of extending the amplitude two-form A to the three 

( , \ 



component amplitude vector A and is thus referred to as the extension matrix, s 



and p = 
given by 

and 



Py 



are unit vectors that embody the vectorial nature of the beam and are 
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s = e u 



kp k x 



(7) 



(8) 



k k 

respectively. Because s, p, and k are mutually orthogonal, Eq. together with Eqs. 
(SJ)-^SD, constitutes an integral transformation solution to the wave equations (CQ) and (J2j) . 
It deserves mentioning that the amplitude two-form A of the angular spectrum can be 
arbitrarily chosen. 

Eqs. (j7]) and (jSJ) show that only element p x of the extension matrix produces the lon- 

A 



gitudinal component. If A 







we arrive at the transverse-electric beam mode 







The principle of duality predicts that A = ( | corresponds to the transverse-magnetic 
beam mode. 

In this paper, we consider only the following amplitude two-form that is independent of 
the azimuthal angle (p, 



A 



A{k p ) = lA{k p ), 



(9) 



where l s and l p are constants, I 



describes the polarization state of the angular 



spectrum and is assumed to satisfy the normalization condition \l s \ 2 + \ l p \ 2 = 1, and A(k p ) is 
the amplitude distribution of the angular spectrum. Let us consider the following Gaussian- 
like distribution function, 

A{k p ) = A ex V (-y^f) A m (K), (10) 

where A Q is a constant and Wq is the characteristic width in the transverse dimension. The 
Gaussian factor guarantees that the beam carries finite energy. Different choices of the 
modulation factor, A m , will correspond to different kinds of beam as will be shown below. 
For the sake of simplicity, we discuss only the paraxial beam, for which condition 

A0=-*-<l (11) 

kw 



holds 



161 ]. where A9 is half the divergence angle that is determined by the Gaussian factor 



in Eq. fflOl) . Under this paraxial condition, Eq. ([3]) can be rewritten as 

E(x) = — / k p dk p \ dipA(k p ,ip)exp(ik-x), (12) 

^ J-oo Jo 

where the integration limits have been extended to ±oo for the variable k p . And the Gaussian 
factor in Eq. (flOl) indicates that the quantity % in the extension matrix can be regarded 
as a small number in comparison with unity when integral ffl2|) is considered. This explains 
why the longitudinal component of a Gaussian-like paraxial beam is of the first order in 
comparison with the zeroth-order transverse component [if]]. Substituting Eqs. (T4|) and 
(|6"j)- (lTUj) into Eq. (fT2"|) and with the help of the following expansion, 

oo 

exp(ip cos ip) = i m J m {p) exp(im^), (13) 

m=— oo 

where J m 's are the Bessel functions of the first kind, we obtain for the electric-field vector, 
E(x) = [i(l s e^ + l p e r )E T (r,x) - l p e x E L (r, x)] exp(ikx) , (14) 



where 



/•oo 

E T (r,x)= A'ik^Mrk^kpdk,, (15) 
Jo 

f°° k 

E L (r,x) = J^ -^A'(k p )J (rk p )k p dk pj (16) 
A\k p ) = A exp (~Y k2 p) A ™( k p)' 



w = Wr. 1 + i — 

and xr = kwQ, which represents the diffraction length. In deriving Eq. flHj) , we have also 

k 2 

made (i) the paraxial approximation [17| k x ~ k — in the exponential factor exp(zk • x), 



2fe 



(ii) and the zeroth-order approximation p p = -^e p « e p in the extension matrix. 

Eq. (j!4p describes bound beams that are axisymmetric with respect to the propagation 
axis not only in the polarization but also in the complex amplitude. By "axisymmetric" we 
mean "invariant" under arbitrary rotation about the axis. The first term on the right side 
is the transverse component. Its amplitude, given by Eq. (fT5l) . is the Hankel transformation 
[l8| of order one of the function A' and is of the zeroth order. The second term is the 
longitudinal component. Its amplitude, given by Eq. (Tl6l) . is the Hankel transformation of 
order zero of the function %A' and is therefore of the first order, ~ ^. So the longitudinal 



component is much smaller than the transverse component 16]. Neglecting the small lon- 
gitudinal component, the beam is dark on the axis r = and is locally polarized elliptically 
with the same polarization state as that of the angular spectrum, I. 

Let us now look at a few examples, paying our attention mainly to the amplitude of the 
transverse component. 

1 Doughnut modified- Bessel- Gaussian vector beams I For the simplest modification factor, 

A m = 1) 

we obtain for the amplitude of the transverse component, 



/27T / r 2 

E T (r,x) = — — -A rexp 



Aw 3 V Aw 2 



^>2 \ / ^,2 



Aw 2 J " \Aw' z 

where Iq and I\ are the modified Bessel functions of the first kind. Due to the linear factor 
r, the beam is dark on the axis. In addition, there is only one bright ring in the transverse 
intensity distribution. So this is a doughnut beam. 

2 Doughnut modified-Bessel- Gaussian vector beams II If we choose for the modification 
factor, 



A m = ex P [-^{3 2 ) Ji(/3w 2 k p ), 
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where (3 is a constant, we have 



TP ( \ A ° 

E T {r,x) = -^-exp 



w- 



cxp 



1+4 



X 



2w< 



This is also a doughnut beam. The radius of the bright ring expands with the increase of 
the value of (3. But the width of the ring changes little. In Fig. [T] is shown the dependence 
of the transverse intensity, I = |-&r| 2 ; on the radial coordinate r at the focal plane x — 0, 
where kwo = 1000, the intensity is normalized to unity, r is in units of wavelength A, the 
solid curve is for [3 = 0.002fc, and the dashed curve is for (3 = 0.005A;. 
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FIG. 1: Dependence of normalized transverse intensity on the radial coordinate at the focal plane, 
where kwo = 1000 and r is in units of wavelength A. Solid: = 0.002&; dashed: j3 = 0.005/c. 



To the best of our know, this is the first time to observe theoretically the above-mentioned 
2 kinds of vector beam. Since the modified-Bessel-Gaussian scalar beams show an elongated 
diction-free region Q in —son with the fnnda m enta. Gans,an bea,n, the p r op- 
agation properties of the modified-Bessel-Gaussian vector beams deserve investigation in 
detail. This is beyond the scope of this paper and will be presented elsewhere. 

3 Bessel-Gaussian vector beams With a modification factor containing the modified Bessel 
function of the first kind of order one, 



A 7 



exp ( ) h((3w' k p ), 



we find 



A 



E T (r,x) = ^J 1 (^/3r)exp 



w- 
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If I 




, we arrive at the azimuthally polarized J\ Bessel- Gaussian vector beam, the 



same as was obtained by solving the paraxial wave equation 9j. 

4 Laguerre- Gaussian vector beams Furthermore, with a modification factor containing 
the associated Laguerre polynomial, 



where a is a constant, we have 

E T (r,x) -- 



A m 


_ k p j i 


A> 
kw A 


/ a 2 
\ w z 


exp I 


' r 2 \ 


2w 2 ) 



r x 



Li 



a 



w 2 + a 2 2w 2 



This kind of vector beam includes the axisymmetric Laguerre-Gaussian vector beam dis- 
cussed in Ref. [12 ]. 



When the amplitude two-form A depends on the azimuthal variable in the wavevector 
space, the vector beam produced by Eq. ([3]) will in general be no longer axisymmetric. 
Suitable choices will lead to cylindrical vector beams of topological charges 
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12 



20j. 
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